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Abstract 



We employ the familiar canonical quantization procedure in a given cosmological setting to 
argue that it is equivalent to and results in the same physical picture if one considers the deforma- 
tion of the phase-space instead. To show this we use a Probabilistic Evolutionary Process (PEP) 
to make the solutions of these different approaches comparable. Specific model theories are used 
to show that the independent solutions of the resulting Wheeler-DeWitt equation are equivalent 
to solutions of the deformation method with different signs for the deformation parameter. We 
QQ I also argued that since the Wheeler-DeWitt equation is a direct consequence of diffeomorphism 

■ invariance, this equivalence is only true provided that the deformation of phase-space does not 

[ break such an invariance. 
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1 Introduction 

X 

^ • Standard cosmological models based on classical general relativity have no convincing and precise 

^ answer for the presence of the so-called "Big-Bang" singularity. Any hope of dealing with such 

singularities would be in vein unless a reliable quantum theory of gravity can be constructed. In the 
absence of a full theory of quantum gravity, it would be useful to describe the quantum states of the 
universe within the context of quantum cosmology, introduced in the works of DeWitt [1] and later 
Misner [2]. In this formalism which is based on the canonical quantization procedure, one first freezes 
a large number of degrees of freedom and then quantizes the remaining ones. The quantum state of 
the universe so obtained is then described by a wave function in the mini-superspace, a function of 
the 3-geometry of the model and matter fields presented in the theory, satisfying the Wheeler-DeWitt 
(WD) equation. In more recent times, the research in this area has been quite active with different 
approaches [3], see [4] for a review. Interesting applications can be found in [5] where canonical 
quantization is applied to many models with different matter fields as the source of gravity. 

An important ingredient in any model theory related to the quantization of a cosmological model 
is the choice of the quantization procedure used to quantized the system. As mentioned above, 
the most widely used method has traditionally been the canonical quantization method based on 
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the WD equation which is nothing but the apphcation of the Hamiltonian constraint to the wave 
function of the universe. However, one may solve the constraint before using it in the theory and in 
particular, before quantizing the system. If we do so, we are led to a Schrodinger type equation where 
a time reparamctcrization in terms of various dynamical variables can be done before quantization 
[6, 7, 8]. A particularly interesting but rarely used approach to study quantum effects is to introduce 
a deformation in the phase-space of the system. It is believed that such a deformation of phase-space 
is an equivalent path to quantization, in particular to canonical and path integral quantizations [9]. 

An important question then naturally arises in applying the various quantization methods to 
cosmological models, that is, if they are equivalent. To at least partially answer the above question, 
we propose to quantize some simple cosmological models, namely the de Sitter, dusty FRW, FRW 
with radiative matter and Bianchi type I, using the methods described above. First we introduce a 
deformation to the Poisson algebra of the corresponding phase-space of the models. This will lead 
us to a deformed Hamiltonian from which the equations of motion can be constructed. We will show 
that the presence of the deformation parameter in the solutions can be interpreted as a quantum 
effect. This is done by comparing the resulting solutions with that of the other quantization method 
which is nothing but the usual canonical WD approach. Indeed, we will show that both of these 
quantization methods have the same physical interpretation for our chosen models. 

A remark on the WD approach is that, as is well known, in canonical quantization the evolution of 
states disappear since in all diffeomorphism invariant theories the Hamiltonian becomes a constraint. 
Indeed, the wave function in the WD equation is independent of time, i.e. the universe has a static 
picture in this scenario. On the other hand in the deformed phase-space method, time does appear 
and so the evolution becomes meaningful. Therefore, the problem of the evolution of states which 
is a major problem in quantum cosmology may also be addressed in this approach. To propose a 
possible solution to this problem we will use what we call the "Probabilistic Evolutionary Process" 
which was introduced in [10, 11]. This mechanism is briefly studied in the next section because of 
its crucial role in this work. In section 3 we review the quantization of the de Sitter and dusty FRW 
cosmological models inspired by the Deformed Special Relativity (DSR), discussed in [12, 13]. The 
canonical quantization for both models will be studied in section 4. We will close the paper with a 
discussion and comparison of the results. 

2 Probabilistic Evolutionary Process (PEP) 

To quantize a classical model the following procedure is usually followed. The classical Hamiltonian 
is written in its corresponding operator form and the resulting Schrodinger equation obtained after 
quantization, i.e. i'h^^ = H^, becomes the relevant equation to describe time evolution of the 
quantum states. However, in diffeomorphism invariant models the Hamiltonian becomes a constraint, 
H = 0, and therefore does not provide for the evolution of the correspond states. This means that 
in such models all the states are stationary. One of the common examples of this situation is general 
relativity which is employed to investigate the evolution of the cosmos. In quantum cosmology, the 
Schrodinger equation becomes the Wheeler-DeWitt equation, H'^ = 0. As is well known, quantum 
cosmology suffers from a number of problems, namely the construction of the Hilbert space for defining 
a positive definite inner product of the solutions of the WD equation, the operator ordering problem 
and most importantly, the problem of time; the wave function in the WD equation is independent 
of time, i.e. the universe has a static picture in this scenario. This problem was first addressed in 
[1] by DeWitt himself. However, he argued that the problem of time should not be considered as a 
hinderance in the sense that the theory itself must include a suitable well-defined time in terms of its 
geometry or matter fields. In this scheme time is identified with one of the characters of the geometry, 
usually the scale factor and is referred to as the intrinsic time, or with the momentum conjugate to 
the scale factor or even with a scalar character of the matter fields coupled to gravity in any specific 
model, known as the extrinsic time. 

In general, the crucial problem in canonical quantum gravity is the presence of constraints in the 
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gravitational field equations. Identification of time with one of the dynamical variables depends on 
the method we use to deal with theses constraints. Different approaches arising from these methods 
have been investigated in details in [6]. As discussed in [6], time may be identified before or after 
quantization has been done. There are approaches, on the other hand, in which time has no funda- 
mental role. The problem is how one can describe the evolution of the universe since observations 
show that the universe is not presently in a stationary state. In a previous work [10, 11], we intro- 
duced a mechanism which we have called the Probabilistic Evolutionary Process (PEP), based on 
the probabilistic structure of quantum systems, to provide a sense of the evolution embedded in the 
wave function of the universe. This is based on the fact that in quantum systems the square of a 
state defines the probability, Va = |^'(a)p. The mechanism introduced as PEP says that the state 
*a makes a transition to the state '^a+da if their distance, c/a, is infinitesimal and continuous, and 
that the higher the value of the transition probability^ the larger the value of Va+da — 'Pa^ ■ The 
mechanism for transition from one state to another is through a small external perturbation^. To 
make the discussion above more clear, we take an example which we shall encounter later on but will 
present the result in the form of the following figure. In what follows, we shall focus on the proba- 
bilistic description of our quantum states provided by figure 1 without worrying about the details of 
the model of which the figure is a result. 




Figure 1: Figure used to explain the idea of PEP. 

Let us take a specific initial condition, say a = 2.5, corresponding to the point P. Then, PEP 
states that the system (here specified by the scale factor a) moves continuously to a state with higher 
probability and thus P moves to the right to reach the point Q, a local maximum. Here, since Q 
locally has the maximum probability the system stays at Q. This means that the scale factor becomes 

A PEP 

constant as the time passes . We show this transition by P — >Q. Now let the initial condition be 

pj^p PEP 

the point R. Then we have R — >Q, and so on. Note that R — >S is possible but it has much 

smaller probability compared to the transitions R^-^Q. We note that the transitions R^^S and 
PEP 

S — > T may be interpreted as tunnelling precesses in ordinary quantum mechanics in the sense that 
the probability of being at S is zero. It means that PEP can reproduce tunnelling processes but with 
a very small probability. 

In the following we will use PEP to describe the evolution of quantum cosmological states in the 

^This transition probability can play the role of the speed of transition i.e. the higher the probability of transition 
the larger the speed of transition. 

^ Since there is no constraint on the positivity of Va+da — 'Pa then PEP can describe a tunneling process too. 

^Certainly, in quantum cosmology, the universe is considered as one whole [14] and the introduction of an external 
force is irrelevant. However, because of the lack of a full theory to describe the universe, these small external forces are 
merely used to afford a better understanding of the discussions presented here. 

''a perturbation around the local ma:ximum is acceptable as described before. 
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canonical method of quantization. We insist that PEP plays a crucial role in the interpretation of 
the states in canonical quantization and allows them to be compared to the states resulting from 
quantization by deformation of the phase-space structure. 



3 Phase-space deformation: a procedure for quantization 

It has long been argued that a deformation in phase-space can be seen as an alternative path to quan- 
tization, based on Wigner quasi-distribution function and Weyl correspondence between quantum- 
mechanical operators in Hilbert space and ordinary c-number functions in phase-space, see for example 
[9] and the references therein. The deformation in the usual phase-space structure is introduced by 
Moyal brackets which are based on the Moyal product [15]. However, to introduce such deformations 
it is more convenient to work with Poisson brackets rather than Moyal brackets. 

From a cosmological point of view, models arc built in a minisuper-(phase)-space. It is therefore 
safe to say that studying such a space in the presence of deformations mentioned above can be 
interpreted as studying the quantum effects on cosmological solutions. One should note that in 
gravity the effects of quantization are woven into the existence of a fundamental length [16]. The 
question then arises as to what form of deformations in phase-space is appropriate for studying 
quantum effects in a cosmological model? Studies in noncommutativc geometry [17] and generalized 
uncertainty principle (GUP) [18] have been a source of inspiration for those who have been seeking 
an answer to the above question. More precisely, introduction of modifications in the structure of 
geometry in the way of noncommutativity has become the basis from which similar modifications 
in the phase-space have been inspired. In this approach, the fields and their conjugate momenta 
play the role of coordinate basis in noncommutativc geometry [19]. In doing so an effective model 
is constructed whose validity will depend on its power of prediction. For example, if in a model 
field theory the fields are taken as noncommutativc, as has been done in [19], the resulting effective 
theory predicts the same Lorentz violation as a field theory in which the coordinates are considered 
as noncommutativc [20]. As a further example, it is well known that string theory can be used to 
suggest a modification in the bracket structure of coordinates, also known as GUP [18] which is used 
to modify the phase-space structure [21]. Over the years, a large number of works on noncommutative 
fields [15] have been inspired by noncommutative geometry model theories [17]. 

In this paper we study the effects of the existence of a fundamental length in a cosmological 
scenario by constructing a model based on the noncommutative structure of the Deformed (Doubly) 
Special Relativity [22] which is related to what is known as the K-deformation [23]. This way of 
introducing noncommutativity is interesting because of its compatibility with Lorentz symmetry, as 
is commonly believed [23, 24]. The K-deformation is introduced and studied in [25]. The «;-Minkowski 
space [26] arises naturally from the K-Poincare algebra [22] such that the ordinary Poisson brackets 
between the coordinates are replaced by 



where k is the deformation (noncommutative) parameter which has the dimension of mass k = ei~^ 
when c= h = 1, and e = ±1 [27] such that k and i arc interpreted as dimensional parameters which 
are identified with the fundamental energy and length, respectively. As mentioned above, one can 
change the structure of the phase-space based on equation (1). Here we will examine a new kind of 
modification in the phase-space structure inspired by relation (1), much the same as has been done 
in [15, 21, 28]. In what follows we introduce noncommutativity based on K-Minkowskian space and 
study its consequences on the de Sitter and dusty FRW cosmologies. 

Let us start by briefly studying the ordinary, spatially flat FRW model where the metric is given 



(1) 



K 



by 



ds^ = - 



N'^{t)dt^ + a^{t){dx^ + dy^ + dz'^), 



(2) 
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with N(t) being the lapse function. The Einstein-Hilbert Lagrangian with a general energy density 
V{a) becomes 

jC = ^{R[g]-V{a)) 

= -m^^aa'^ - Na^V{a), (3) 

where R[g] is the Ricci scalar and in the second line the total derivative term has been ignored. The 
corresponding Hamiltonian up to a sign becomes 

^0 = ^Na-Va-Na^V{a). (4) 

Here, we note that since the momentum conjugate to N(t), tt = |^ vanishes, the term Att must be 
added as a constraint to Hamiltonian (4). The Dirac Hamiltonian then becomes 

n = ^Na-^pl- Na^V{a) + Xtt. (5) 

To introduce noncommutativity one can start with 

{N'{t),a'{t)} = £a'{t). (6) 

This is similar to equation (1) since one can interpret N{t) and a{t), appearing as the coefficients of 
dt and dx, in the same manner as xq and Xi appearing in (1) respectively. For this reason we shall 
call it the K-Minkowskian-minisupcr-phasc-space. In this case the Hamiltonian becomes 

n'o = ^N'a'-'p'l-N'a"v{a'), (7) 

where the ordinary Poisson brackets are satisfied except in (6). To move along, one introduces the 
following variables [12, 13, 29] 

' N'{t) = N{t)-£a{t)pa{t), 
a'it) = ait), (8) 

. P'ait)=Pa{t). 

It can be easily checked that the above variables will satisfy (6) if the unprimed variables satisfy 
the ordinary Poisson brackets. The term —ia{t)pa{t) may be looked upon as a direct consequence 
of a phase-space deformation of relation (6) which, as has been suggested, could originate from 
string theory, noncommutativc geometry and so on, see [9, 28, 30]. With the above transformations, 
Hamiltonian (7) takes the form 

= 2i^«~ VI - Na'Via) - ^£pl + £a^V{a)pa. (9) 

It is clear that the momentum conjugate to N{t) does not appear in (9), i.e. tt = is a primary 
constraint. It can be checked by using Legendre transformations that the conjugate momentum 
corresponding to N{t), vr = J^, vanishes. It is therefore necessary to add the term Att to Hamiltonian 
(9) to obtain the Dirac Hamiltonian 

= ^^""''pI - Na^V{a) - j^ipl + £aV(a)p„ + Att. (10) 
The equations of motion resulting from Hamiltonian (10) are 

d = {a,n'"'} = ^Na-'pa-^epl + ea''V{a), 

Pa = {Pa, H""} = ^^a'Va + SNa^V{a) + Na^V'{a) - 4^aV(a)p„ - ^a^F'(a)p„, 

AT = {N,n''^} = X, 

n = {n,nn = -^a-Va + a'V{a), (11) 
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Figure 3: The scale factor for t = —0.0001 in the de Sitter model. The initial conditions are Ci = -0.0001 and A = 6 
for the left graph and Ci = 1 for the right one. 



where a prime denotes differentiation with respect to the argument. The requirement that the primary 
constraints should hold during the evolution of the system means that vr = {tt,^"'^} = 0. If is now 
calculated from the secondary constraint, tt = 0, and the result is substituted in the first equation in 
(11) one obtains 

d + 2ia'V{a) = N^^a'^V{a). (12) 
It is easy to check that the above equation is consistent with the second equation in (11) as well. 



3.1 The de Sitter model 

As is well known, in the de Sitter model V{a) = 2A and equation (12) becomes 



d + UAa"^ = N^^Aa^, (13) 

which is compatible with other equations. The solution of the above equation, the scale factor, for a 
constant N can be written as^ 

1 / ATeV^^i \^ 
a(t) = 36 ^ ^ , (14) 

where Ci is a constant of integration. It is seen that the usual de Sitter solution (without any 
deformations) is recovered in the limit ^ ^ 0. The scale factor calculated in (14) and its behavior are 



^This means that wc restrict ourselves to a certain class of gauges, namely N = const., which is equivalent to the 
choice A = in the equation of motion (11). 
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shown in figure 2 for ^ > and in figure 3 for ^ < 0. In what foUows, we discuss the different figures 
separately. This may seem irrelevant here since the behavior of the scale factors for different choices 
are similar to each other. However, in the next case when we study the dust model, the solutions are 

essentially different, demanding separate discussions on their behavior. The same is also true when 
we study the relation of the deformed phase-space solutions to that of the canonical method at the 
end of next section. 

The left plot in figure 2 shows that the scale factor has started from an initial value and becomes 

constant in the end. During this passage the behavior of the scale factor is monotonically increasing. 
The right plot in figure 2 on the other hand, shows a conversing behavior and predicts a monotonically 
decreasing scale factor. Note that in this case the scale factor becomes constant in its final state too. 
The scale factor represented by the left plot in figure 3 has the same behavior as that in the left 
plot in figure 2. In figure 3 one must restrict the discussion to the second part of the right plot only 
since the denominator of the scale factor in (14) cannot be zero for physical reasons. This means, for 
example, that the universe has come into being at t = 1.5 in our case. Now the behavior of the scale 
factor is the same as that of the scale factor in the right plot of figure 2. Note that in all the figures 
we have taken A'^ = 1 which only changes the values of the initial and final scale factor without any 
change in the behavior''. It should be mentioned here that these different initial and final values are 
crucial to our discussion on the equivalence between these two different approaches. 



3.2 The dust model 

In this case, substituting V{a) = poa~^ in equation (12) results in 



a + 2ipQa = 



(15) 



The solution for N = const, is given by 

a{t) 



2 



for ^ > and 



a{t) 



2x33 



2x33 



1 - C2e-^^''o* 



(16) 



(17) 



for i < 0, where C2 is an integration constant which, for the second solution (17), must satisfy 
C2e3|^lwt > 1^ indicating that C2 > 1. Fi gures 4 and 5 show the behavior of the scale factors for 
£ > and £ < respectively. The scale factor presented by the left plot in figure 4 suggests that 
the universe starts from an initial state and reaches a constant final state due to a monotonically 
increasing behavior. As can be seen from the right plot in figure 4, the converse is also true. It is 
worth to consider that only for C2 = 1 the limitation i results in exactly the same behavior due 
to its counterpart in absence of any deformations^. 

As mentioned in the previous section, the FRW dust model represents a completely different 
behavior for ^ < 0, which is represented in figure 5. In this case the scale factor monotonically 
increases without reaching a final constant value. 



This is compatible with the spirit of gauge invariance, that is, there is no physical difference between multiple gauge 
fixing. The measured value of the system parameters could be different for different gauges but as mentioned before, 
they have the same behavior. 

'^It does not make other choices for C2 irrelevant since the limiting process is yet a questionable matter. Specifically, 
the solution which is represented in the right figure in figure 4, has not any correspondents in non-deformed phase-space 
and it is a prediction of quantum cosmology only. 
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Figure 4: The scale factor for I = 0.01 in the dust model. The initial conditions are C2 = 1 and po = 1 for the left 
graph and C2 = — 1 for the right one. 

a(t) 




Figure 5: The scale factor for I = —0.01 in the dust model. The initial conditions are C2 = 1 and A = 6. Note that 
in this case, for all the valid values of C2 {C2 > 1) we have the same behaving graphs. 



3.3 The radiation model 

This case is more interesting since it aims to describe the early phases of the Universe, when quantum 
geometrical properties are expected to dominate. Mathematically, for radiation era V{a) = poO'~^ 
which makes equation (12) 



d + 2lpQ 

The solution of the above equation reads as 



6 



a{t) 



i + w 



V6N 



J 



(18) 



(19) 



for £ > and it reads 



a{t) 



V6N 



1 + w - 



(20) 



and for ^ < where C3 is an integration constant, W{z) is Lambert W function and gives the principal 
solution for w in z = we!" . Figures 6 and 7 show the behavior of scale factor for ^ > and £ < 0, 
respectively. 



4 Canonical quantization: the WD equation 

We now focus attention on the study of the quantum cosmology of the models described above. For 
this purpose we quantize the dynamical variables with the use of the WD equation, that is, = 0, 
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Figure 6: The scale factor for I — 1, po = 1 and C3 = 1 in the radiation model is shown in the left figure. In the right 
figure the scale factor for C3 = — 1 is shown. 
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Figure 7: The scale factor for £ — —1, po = 1 and C3 = 1 in the radiation model is shown in the left figure. In the 
right figure the scale factor for C3 = 1 is shown. Note that since in the metric (2), a^{t) is appeared so the absolute 
value of the scale factor (20) is meaningful. 

where Ti is the operator form of the Hamiltonian given by equation (4), and ^ is the wave function 
of the universe, a function of the scale factor and the matter fields, if they exist. 



4.1 The de Sitter model 

Let us set V{a) = 2 A. Then the corresponding WD equation becomes 



1 -1 2 

— a Vn 

24 



2Aa^ 



^{a) = 0. 



(21) 



Choice of the ordering a = paa ^pa to make the Hamiltonian Hermitian and use of [a,pa\ = ih 
and Pa = —ifida results in 



with solutions 



/4aVA\ 



(22) 



(23) 



where J is the Bessel function and ci and C2 are arbitrary constants. The constants must be so chosen 
as to satisfy the initial conditions. However, in what follows, our discussions are independent of the 
values of these constants. Here we show and emphasize the different behavior of both solutions. To 
do this we need to calculate the probability density |^'(a)p. In figures 8 and 9 the probability density 
for C2 = and ci = are shown respectively. 

Now, we start form a chosen initial condition, say, a = 1. The PEP procedure may now be 
employed to describe the physical interpretation of the resulting states. In figure 8, let the initial 
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Figure 8: The probability density for the Ji term with A = 0.0001. 
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Figure 9: The probability density for the J_ i term with A = 0.0001. 



state be at point P. One then expects, keeping in mind the discussion presented previously, the 

PEP 

transition P — '>Q to appear, with the fohowing physical interpretation. The point P shows an initial 
state and the point Q shows the final state. The universe begins its evolution by a monotonically 
increasing scale factor to reach the point Q. Point Q is a state with a greater scale factor and is finite. 
Note that since the point Q is a local maximum, PEP predicts Q — > Q which means that the scale 

PEP 

factor remains in this final state. The transition P — >Q then means starting from an initial state, 
increasing monotonically, and finally arriving at a constant final state. This behavior is completely 
similar to the behavior of the scale factor in the previous section, presented in the left plot in figure 
2. 

Since for the other cases the details are similar, in the following we mention only the correspon- 

PEP 

dence between the solutions without any detail. In figure 8, the transition R — >Q is similar to the 

PEP 

right plot in figure 2. On the other hand in the same figure, transition T — )-U is similar to that of the 

left plot in figure 2 with a difi^erent initial condition, that is, with an appropriate Ci. The transition 
PEP 

Q — >Q represents the scale factor (14) with Ci = 0. In figure 9, similar discussions as above apply 
and the transitions Q^^P and T^-^S have the behavior and interpretation compatible with the 

PEP 

right plot in figure 3 with different initial conditions. Again, R — >S is similar to the left plot in 
figure 3. 

PEP 

Now, an interesting question arises as to the meaning of P — >P in this case. One possible answer 
would be that only superpositions of the solutions are appropriate for comparison with the previous 
section results, e.g. ci = C2 = 1 and ci = — C2 = 1. 

It is appropriate at this point to discuss different values of the final scale factor predicted by 
the WD equation. Obviously, since the maximum probabilities in figures 8 and 9 occur at different 
values then the final scale factor assumes different values too. This prediction is compatible with 
the deformed phase-space method due to the appearance of the lapse function, A*", as a constant 
parameter in relation (14) controlling the behavior of the scale factor. It is worth mentioning that 



10 



Probability 



Probability 




V U V u 



Scale factor 




0.5 1 1.5 



Scale factor 



Figure 10: The probability density for the Bi' term with po = 1. 
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Figure 11: The probability density for the Ai' term with po = 1. 

in both approaches there are two parameters that control the behavior of the scale factor; Ci in the 
deformed phase-space method and also A'^ which is fixed by the maximum value of the scale factor 
and its initial values in each region of the WD solution. Note that each region is an interval between 
two minima that contains one maximum. For example in figure 6 the region containing points P and 
Q is one region and that containing T and U is another region and the same is true for other figures. 



4.2 The dust model 

In this case we set V{a) = poa~^ in Hamiltonian (4). The corresponding WD equation becomes 

1 



With the above ordering, the WD equation is 



■^{a) = 0. 



(24) 



dl^{a) - a-^da^ia) + ^a*(a) = 0, 



with solution 



*(a) = c'lAi' 



+ c'nBi' 



(25) 



(26) 



where c'^ and are integration constants and Ai' and Bi' are derivatives of the Airy functions Ai 
and Bi with respect to a respectively. In this case the corresponding probability densities |*(a)p are 

shown in figures 10 and 11 where c'l = and c'2 = respectively. In the dust model two different 
solutions show completely different behaviors. One, figure 10, has local maxima but the other, figure 



.PEP, 



11, is monotonically increasing. In figure 10, the transition Q — >R is compatible with the left plot 

PEP 

in figure 4. Transition S — >R however, is compatible with the right plot in figure 4 with transitions 
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Figure 12: The probability density for the J-Bessel term with po = 1 and h = 1. 
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Figure 13: The probability density for the F-Bessel term with po = 1 and h = 1. 



Scale factor 



P^-^P and R^-^R being compatible with scale factor (16) with different initial conditions and 
constants. 

PEP 

It is worth noting that the transition P — >Q in figure 11 means that no matter what the initial 
state is, the scale factor moves to a larger value and continues to grow since in this case no local 
maximum exists at all. This behavior is completely in agreement with the behavior of the scale factor 
represented in figure 5. Note that a similar discussion like that of the last paragraph of the previous 
subsection goes for the dust case. Note that for figure 11 there is only one region. 

4.3 The radiation model 

In the radiation case the Hamiltonian (4) with V{a) = /Oo^"^ results in the following WD equation 

1 



24 



^-(0) = 0. 



By using the mentioned ordering it reduces to a differential equation as follow 

24po 



with solution 



d^^ia) - a-'da'i^ia) + — ^a^'(a) = 0, 



^{a) = c[aJi {2VQha) + c^aYx (2V6ha 



(27) 



(28) 



(29) 



where c" and C2 are integration constants and J and Y are Bessel functions. The behavior of the 
above functions are represented in figures 12 and 13. The discussions on comparison between quantum 
cosmological solutions and their correspondence from deformed phase-space formalism, i.e. figure 6, 
are the same as previous models, cosmological constant and dust models. To prevent making the 
paper boring we omit discussing for this case. 
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5 Bianchi type I Model 



In this section we consider a more complicated model i.e. Bianchi model with different matter fields 
and compare its quantum behavior that is modeled by a kind of phase-space deformation and also 
canonical quantization. We separate Bianchi models from previous ones since there is more than one 
variable exists that makes our conjecture more reasonable. 

5.1 Phase-space deformation 

Let us consider a cosmological model in which the spacetime is assumed to be of Bianchi type I whose 
metric can be written as 

= _iv2(t)dt2 + e^^We^/^^^Wdx^dx^ (30) 

where N(t) is the lapse function, e"*^*-* is the scale factor of the universe and fiij{t) determine the 
anisotropic parameters v{t) and w{t) as follows 

Pij = diag [y + \/^w, v - ^/^w, -2v^ . (31) 

To simplify the model we take u; = 0, where is equivalent with a universe with two scale factors in 
the form 

ds^ = -N'^{t)dt^ + a^{t){dx^ + dy"^) + c^{t)dz^. (32) 

The anisotropy in the above metric can achieved by introducing a large scale homogeneous magnetic 
field in a flat FRW spacetime. Such a magnetic field results in a preferred direction in space along the 
direction of the field. If we introduce a magnetic field which has only a z component, the resulting 
metric can be written in the form (32) where there are equal scale factors in the transverse directions 
X and y and a different one, c(t), in the longitudinal direction z. The Hamiltonian for gravity coupled 
to a perfect fiuid with equation of state p = 7p is 

-H = (-pI + pI) + NM^e-^-^- + Att, (33) 

where, My is a model dependent constant. To introduce noncommutativity one can start with 

{Ar'(t), e"'+"'} = ^e"'+"', {N'{t), e"'-^'''} = ^e"'-^"', (34) 
where are the two dimensional generalization of relation (6). The Hamiltonian of this model becomes 

■H'o = 2i^'e-3-' (-p'l+p'l) + N'M^e-'-ru'^ (35) 
Now we introduce as before 

' N'{t) = N{t) - aip.^{t) - /3epy{t), 

< u'{t) = u{t), v'{t) = v{t), (36) 
. Pu{t)=Pu{t), p'^{t)=py{t) 

where a + /3 = 1. The deformed Hamiltonian by imposing the above transformations takes the form 

H'"'=[Nit)-a£pu{t)-P£p^{t)] 



-37M 



+ Att. 



(37) 
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Figure 14: Left plot is shown for j = —1 (de Sitter) and right one for 7 = (dust). We take a = l3 = 1/2, pov = 1, 
A'^ = 1, ^ = 1, Mj = 1 and initial conditions u{t = 0) = v{t = 0) = 1, Pu{t = 0) = 0. The dashed lines represent 
M-behavior and the solid line u-behavior in both figures. 



The classical dynamics is governed by the Hamiltonian equations, that is 



u = {n,-?^"^} = -^e-3>„ [N{t) - a£pu{t) - f3£p,{t)] - ai 



J_ -3m 

24^ 



J_ -3« 
12^ 



p, [N{t)-alpu{t)-pip.o{t)]-pi 



Lg-3M 



24 



Pu = {Pu^-W] = [N{t) - aipuit) - Plp.it)] 
Pv = {Pv,n'''^} = 0, 



1 -3u 
8^ 



-3'yu 



(38) 



N = {N,n''^} = A, 

_ ^ = {vr, n^-} = ^e-3« {-pI + pI) + M^e-3T-. 

The requirement that the primary constraints should hold during the evolution of the system means 
that 77 = {vr,^"'^} = 0. Applying this condition into the above equations we get 



u 



J_ „-3u 
12^ 



Pu[N-aipuit)-l3£pov] 



V = 3^e-3>^ [N - atpuit) - ^epov] , 



. Pu = -3(7 - l)M^e-3^« [N - a£pu{t) - pip 



(39) 



Since the analytic solution does not exist for the above equations, their behavior is represented in 
figures 14 and 15 by applying numerical methods. In figure 14 for u, there is a decreasing era before 
steady behavior in contrast to u-behavior. However, the other choices are possible too. For example, 
in figure 15, both of u and v are the same in the existence of a decreasing era in their behaviors. 
These figures are plotted for typical numerical values for the parameters and initial conditions. With 
examining some other initial conditions the behavior is almost repeated. 



5.2 Canonical quantization 

The WD equation for our model reads 



The solutions to the above equation are 



r)2 Pfl 



^'(u, v) = 0. 



M/,(n,z;)=e±*^V±,,/3(4 



7g3MX 



(40) 



(41) 
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u 



1.000 , 
0.999 : 
0.998 : 
0.997 : 
0.996 - 
0.995 - 



Figure 15: This figures are sketched only for dust, 7 = 0. We take a = j3 = 1/2, pov = -1, N ^ 1, £ ^ 1, M-y = 1 
and initial conditions u(t = 0) = v{t = 0) = 1, Puit = 0) = 0. In the left plot it-behavior is represented and w-behavior 
is represented in the right one. 





Figure 16: The right plot is shown |'I'(it, u)p and the left is its corresponding contour plot. It is for de Sitter case 
with = 20. It is taken only two first terms of relation (43) with a constant wight function. 



for 7 



-1 and 



'2M- 



7g3«/2 



(42) 



for 7 = 0, where is a separation constant. We have chosen oscillatory function e^^'^'" since the 
real exponents would lead to exponential increasing wave function for v — >• ±cx) that would not 
show physical behavior. Therefore, we may now write the general solution of the WD equation as a 
superposition of the above eigenfunctions 



+00 



C{v)'^v{u,v)di', 



(43) 



where C{v) can be chosen as a shifted Gaussian weig ht function e""('""^)^ The square of wave 
function for different cases are shown in figures 16 and 17 for de Sitter case and dust respectively. 
In figure 16 point Q is located at a local maximum probability so PEP predicts a transition from 



point P to Q i.e. 



Physically, this transition imposes going u down and going v up and 



after P^-^Q since Q is located at a maximum then Q^-^Q that means the system be steady at Q. 
This behavior is exactly similar to previous result which is represented in the left plot in figure 14. 

PEP 

The similar is true for the left plot in figure 14 and P — >Q transition in figure 17. Since point S 
is located at a local maximum probability so PEP predicts R^-^S then S^-^S. It means there is 
this possibility that both of u and v go down. This later result is completely in agreement with the 
previous subsection results i.e. the represented case in figure 15. 
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Figure 17: The right plot is shown \'^{u, v)\'^ and the left is its corresponding contour plot. It is for dust with M^, — 20. 
It is taken only two first terms of relation (43) with a constant wight function. 

There is a notable difference between Bianchi model and previous ones that is existence of more 
than one direction in configuration space for Bianchi model. This makes some ambiguities in path 
of transition to more probable state in PEP paradigm. But as mentioned roughly in [10] the path 
in these kinds of examples is on the gradient of (hyper-)surface which has the most slope. This 
proposal results in definition of a unique transition path in PEP for more complicated configuration 
spaces. 

6 Discussion 

In this paper we have argued, using familiar examples, that one can draw a completely equiva- 
lent physical interpretation for two different approaches to quantization, namely the usual canonical 
quantization method and the phase-space deformation method, at least within the framework of the 
examples used in this work. A by-product is that when one studies the deformation of phase-space, 
one should realize that this is effectively a quantization procedure and should refrain from using any 
other quantization method simultaneously. However, since in the phase-space deformation method 
for quantization the evolutionary parameter, t, does appear, in contrast to the canonical method, a 
mechanism for retrieving the dynamical information of the system, when canonically quantized, is at 
hand and can be used. This had been done in a previous work [10, 11]. The results show that if the 
equivalence of different methods of quantization discussed above are assumed, then the probabilis- 
tic evolutionary mechanism works consistently and therefore can be used to address the problem of 
evolution in diffeomorphism invariant theories. 

However, such an equivalence cannot hold true in examples where the deformation in phase-space 
is introduced in a Lorentz non-invariant manner. This is not hard to predict since the WD equation 
is a direct result of diffeomorphism invariance and so if a deformation in phase-space breaks such an 
invariance then the results of different quantization methods should be different. In this connection 
we note that K-Minkowski deformation and GUP, for example, preserve Lorentz invariance [24]. 

Another feature which is of interest is that the WD equation is usually a second order differential 
equation and has therefore two independent solutions. However, in the deformed phase-space quantum 
models, the equations of motion are first order, in contrast to the former. The disparity lies in the fact 
that, as discussed above, the parameter of deformation, i, is arbitrary, and hence can be chosen to 
be either positive or negative. It has been shown that different signs of I relate to different solutions 
of the WD equation. 

It is worth mentioning that this paper can be read as a companion for our previous one [10]. In 
[10] we have discussed on the PEP paradigm massively but here to make the current manuscript 
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self-studied we highlight some crucial features of PEP very rapidly. The PEP proposal has some 
roots in the second law of thermodynamics that is a system in an initial state evolves approaching to 
another state in which the final state has a more entropy value or a same one i.e. Sf — Si = AS > 0. 
A vital problem in PEP is the velocity of transition from a state to a more probable state. Actually, 
this question has not any concrete answer yet due to the first steps in understanding PEP. But in a 
model discussed in [10] there is a relation for the velocity of transition what is achieved by comparison 
of PEP's results with Causal Dynamical Triangulation method as a model of quantum gravity. It 
means until now it is an open problem that how the velocity of transition can be addressed by PEP 
itself. 

Finally, it is significant to consider that if a proposal is correct in few examples it does not mean 
its correctness universally and in this paper it is not our claim in any senses. But we have shown that 

there is a kind of correspondence between different methods of quantizations. We suppose, at least, 
this correspondence should be considered until a complete understanding of the universe^. Since in 
lack of a full theory, e.g. quantum gravity, all the approaches are considerable and perhaps each 
approaches sheds some lights on various parts of the universe. This viewpoint results in a proposal 
that a combination of different approaches is a part of full theory^. However, maybe, different 
approaches have an intersection in their territory of validncss that have to be characterized exactly 
to do not result in double consideration, e.g. double calculations. In this paper we have shown that 
two different methods of quantizations can be same as each other at least in some specific examples 
and it must be considered to prevent double quantization. And even more this correspondence can 
shed lights by each methods on the other methods for example it can be interpreted as an evidence 
for correctness of PEP method as a developing in understanding canonical quantization method. 
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